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Abstract 

■^^ I This is a survey of recent results about bipotentials representing multival- 

pLn ■ ued operators. The notion of bipotential is based on an extension of Fenchel's 

(-H I inequality, with several interesting applications related to non associated con- 

stitutive laws in non smooth mechanics, such as Coulomb frictional contact or 
non-associated Driicker-Prager model in plasticity. 

Relations betweeen bipotentials and Fitzpatrick functions are described. 
Selfdual lagrangians, introduced and studied by Ghoussoub, can be seen as 
^S| ■ bipotentials representing maximal monotone operators. We show that bipoten- 

^ I tials can represent some monotone but not maximal operators, as well as non 

monotone operators. 

Further we describe results concerning the construction of a bipotential 
which represents a given non monotone operator, by using convex lagrangian 
I ^ covers or bipotential convex covers. At the end we prove a new reconstruction 

^—v ■ theorem for a bipotential from a convex lagrangian cover, this time using a 

QQ ' convexity notion related to a minimax theorem of Fan. 
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1 Introduction 

In the generalized standard material theory of Halphen and Son [21] any constitutive 
law of a standard material relates a strain rate variable x & X with a stress-like 
variable y S y by using a dissipation potential 0. 

X and Y are topological, locally convex, real vector spaces of dual variables 
X ^ X and y G y, with the duality product (•,•) : X x y ^ R. The dissipation 
potential (/> is a convex and lower semicontinuous function defined on X and the 
associated constitutive law is given by one the following equivalent conditions: 

(a) y G d(f){x), where d(j) is the subdifferential of (j) in the sense of Convex Analysis, 

(b) X G d(j)*{y), where </)* is the Fenchel dual of (j) with respect to the duality 
product, 

(c) (t>{x) + (j)*{y) = {x,y) . 



The constitutive laws of standard materials are also called associated laws. From 
the mathematical viewpoint such laws are cyclically monotone operators. However, 
there are many non-associated constitutive laws, described by a multivalued operator 
T : X ^ 2 which is not cyclically monotone, in some cases not even 
monotone. 

A possible way to study non-associated constitutive laws by using convex anal- 
ysis, proposed first in [33j, consists in constructing a " bipotential" function b of two 
variables, which physically represents the dissipation. See definition 13.11 and the 
section [3] for the introduction into the subject of bipotentials. 

A bipotential function b is bi-convex, satisfies an inequality generalizing Fenchel's 
one, namely ^x £ X,y £ Y, b{x,y) > {x,y), and a relation involving partial subd- 
ifferentials of b with respect to variables x, y. In the case of associated constitutive 
laws the bipotential has the expression b{x, y) = <j){x) + (p* (y) (section [3]) . The graph 
of a bipotential b is simply the set M{b) C X x Y oi those pairs (x, y) such that 
b{x,y) = {x,y). 

A maximal cyclically monotone operator T : X — > 2 is represented by a lower 
semicontinuous and convex "potential" function (p, by a well known theorem of 
Rockafellar. More general, we say that the bipotential b represents the multivalued 
operator T if the graph of T equals M[b). 

There are already many applications of bipotentials to mechanics. Among them 
we cite: Coulomb's friction law [34J, surveyed here in section [9.21 non-associated 
Driicker-Prager [35] (section l9.ip and Cam-Clay models [36] in soil mechanics, cyclic 
plasticity ([34], [3]) and viscoplasticity [24] of metals with non linear kinematical 
hardening rule, Lemaitre's damage law [2], the coaxial laws ( |39j . [44j ) . details in 
sections 18.21 and 19. 3i 

The notion of a bipotential associated to a multivalued operator is interesting 
also from a mathematical point of view. We show that bipotentials are related 
to Fitzpatrick functions associated to a maximally monotone operator. Selfdual 
lagrangians introduced and studied by Ghoussoub jl8j can be seen as bipotentials 
representing maximal monotone operators (section [5|). In section [7] we describe a 
result from [11] which implies that some monotone non maximal operators can be 
represented by bipotentials. 

Other examples of bipotentials come from inequalities. For example, the Cauchy- 
Bunyakovsky-Schwarz inequality can be recast as: if X = y is a Hilbert space then 
the function b : X x X ^^ M. defined by b{x,y) = ||x||||y|| is a bipotential. More 
general, some inequalities involving eigenvalues of real symmetric matrices can be 
put in a similar form, thus providing more non trivial examples of bipotentials. 

In order to better understand the bipotential approach, in [10], [11] we solved two 
key problems: (a) when the graph of a given multivalued operator can be expressed 
as the set of critical points of a bipotentials, and (b) a method of construction of 
a bipotential associated (in the sense of point (a)) to a multivalued, typically non 
monotone, operator. 

At the end of this paper we prove a a new reconstruction theorem for a bipotential 



from a convex lagrangian cover, this time using a convexity notion related to a 
minimax theorem of Fan. 

2 Notations and first definitions 

X and Y are topological, locally convex, real vector spaces of dual variables x G X 
and y G y, with the duality product (•, •) : X x y ^ M. The topologies of the 
spaces X, Y are compatible with the duality product, that is: any continuous linear 
functional on X (resp. 1") has the form x >-^ {x, y), for some y &Y (resp. y ^^ {x, y), 
for some x £ X). 

We use the notation: 1 = M U {+cx)}. 

Given a function (p : X ^ M., the domain dome/) is the set of points x £ X 
with (pix) G M. The polar of (p, or Fenchel conjugate, (p* : y ^ M is defined by: 
(p*{y) = sup{(y,x) - (p{x) \ x G X}. 

We denote by T{X) the class of convex and lower semicontinuous functions (p : 
X ^M. The class of functions (p G T{X) with domcp 7^ is denoted by Tq(X). The 
class of convex and lower semicontinuous i;^ : X ^ M is denoted by T(X, M). 

For any convex and closed set A C X, its indicator function, XA, is defined by 

. . _ f if X G A 

[ +00 otherwise 

Definition 2.1 The graph of an operator T : X ^ 2 is the set G{T): 

G{T) = {{x,y)£XxY :y£T{x)] 
To a graph M C X x Y we associate the multivalued operators: 
X 3 X h^ m{x) = {y G y I (x, y) G M} , 

Y 3 y^ m*{y) = {x G X | (x,y) G M} . 

The domain of the graph M is dom(M) = {x G X | m(x) 7^ 0}. The image of 

the graph M is the set im{M) = {y G y | m*{y) / 0}. 

Definition 2.2 The subdifFerential of a function cp : X ^M. in a point x G dam (p 
is the (possibly empty) set: 

d(p{x) = {ueY \yzeX {z-x,u) < (p{z) - (p{x)} . 

In a similar way is defined the subdijjerential of a function tp : Y ^ M. in a point 
y G domip, as the set: 

dip{y) = {f G X I \/w £ Y {v,w — y) < ip{w) — ip{y)} ■ 

With these notations and definitions we have the Fenchel inequality. 



Theorem 2.3 Let (p : X ^M be a convex lower semicontinuous function. Then: 
(i) for any x £ X,y £Y we have <^(x) + (t)*{y) > {x,y); 
(a) for any {x,y) £ X xY we have the equivalences: 

y G d(t>{x) ^^ X G d(t>*{y) ^^ 0(x) + 0*(t/) = {x,y) . 

Definition 2.4 An operator T : X ^ 2 is monotone if for any {{x,y), {x' ,y') G 
G{T) we have 

{x-x',y -y') > 

A graph M C X x Y is monotone if it is the graph of a monotone operator. The 
graph is maximal monotone (or the associated operator is maximal monotone) if 
for any monotone graph M' C X x Y such that M C M' we have M = M' . 

An operator T : X ^ 2 is cyclically monotone if its graph G(T) is cyclically 
monotone. A graph M is cyclically monotone if for all integer m > and any 
finite family of couples {xj, yj) G M,j = 0,1, . . . ,m, 

m 

(xo - Xrn, y^n) + ^{xk - Xk-i,yk-i) < 0. (2.0.1) 

fc=l 

A cyclically monotone graph M is maximal if it does not admit a strict prolongation 
which is cyclically monotone. 

3 Bipotentials 

Definition 3.1 ^ bipotential is a function b : X x Y ^M with the properties: 

(a) b is convex and lower semicontinuous in each argument; 

(b) for any x £ X,y ^Y we have b{x,y) > {x,y); 

(c) for any {x,y) £ X xY we have the equivalences: 

yedb{-,y){x) ^^ xedb{x,-){y) ^^ b{x,y) = {x,y) . (3.0.1) 

The graph of b is 

M{b) = {{x,y)£XxY I b{x,y) = {x,y)} . (3.0.2) 



An equivalent definition of a bipotential comes out from the following proposi- 
tion. 



Proposition 3.2 A function b : X x Y ^ M is a bipotential if and only if the 
following conditions are satisfied: 

(A) b is convex and lower semicontinuous in each argument and for any x G X, y E 
Y we have b{x,y) > {x,y); 

(Bl) for any y gY , if the function z £ X h^ {^{z, y) — {z, y)) has a minimum, then 
the minimum equals 0; 

(B2) for any x G X , if the function p gY h^ {b{x,p) — {x,p)) has a minimum, then 
the minimum equals 0. 



Proof. Condition (A) is the same as conditions (a),(b) from definition 13. 11 All we 
have to prove is: if the function b satisfies condition (A) then condition (c) from 
definition 13.11 is equivalent with (Bl) and (B2). 

Assume (A) and take x G X, y £ Y such that x £ db{x, •){y). This is equivalent 
with: a; is a minimizer of the function 

z e X h^ {b{z,y) - {z,y)) 

Therefore the statement x G 56(x, ■){y) <^=^ 6(x, y) = {x, y) is equivalent with (Bl). 
In the same way we prove that y G db{-,y){x) -^^ b{x,y) = {x,y) is equivalent 
with (B2). D 

This simple proposition justifies the introduction of strong bipotentials, which 
are particular cases of bipotentials. Conditions (BIS) and (B2S) appeared first time 
as relations (51), (52) pO] . 

Definition 3.3 A function b : X x Y ^ M. is a strong bipotential if it satisfies 
the conditions: 

(a) b is convex and lower semicontinuous in each argument; 

(BIS) for any y gY mi {b{z, y) - {z, y) : z e X} £ {0, +00} ; 

(B2S) for any x £ X inf {6(x,p) — (x,p) .■ p G Y} G {0, +cxd}. 



4 Operators representable by a bipotential 

Definition 4.1 The non empty set M C X x Y is a BB-graph (bi- convex, bi- 
closed) if for all x G dom{M) and for all y G im{M) the sets m{x) and m*{y) are 
convex and closed. 

The following theorem (theorem 3.2 [10]) gives a necessary and sufficient condi- 
tion for the existence of a bipotential associated to a constitutive law M. 



Theorem 4.2 Given a non empty set M G X x Y , there is a bipotential b such 
that M = M{h) if and only if M is a BB-graph. 

The bipotential mentioned in the previous theorem is denoted by 6m and it has 
the expression: 

bM{x,y) = {x,y) +XM{x,y) (4.0.1) 

In the case of a maximal cyclically monotone graph M, by Rockafellar theorem 
( [32] Theorem 24.8.) there is an unique separable bipotential associated to M ( see 
section [6] for the definition of separable bipotentials). With the bipotential given by 
(14.0. ip we have two different bipotentials representing the same graph. Therefore, 
in the larger class made of all bipotentials, in general there is no unicity of the 
bipotential representing a given BB-graph. 

For any BB-graph M, a bipotential b is admissible if M C M{b). Then we 
obviously have b{x,y) < bMix,y) for any (x,y) £ X x Y. In this sense 6m is the 
greatest admissible bipotential for the BB-graph M. 

5 Fitzpatrick functions; selfdual lagrangians 

Let X be a reflexive Banach space and X* its topological dual. The duality product 
between X and X* is the function ir : X x X* -^ M, defined by 7r(x, x*) = {x, x*) = 
x*{x). 

The space X x X* is in duality with itself by the duality product 

{ix,x*),{y,y*)) = {x,y*) + {y,x*) 



Fitzpatrick functions have been introduced in [14]. More on Fitzpatrick functions 
can be found in [1] and in the book [5]. 

Definition 5.1 The Fitzpatrick function associated to a graph M C X x X* is 
the function fM '■ XxX* -^ MU{-|-oo} given by the Fenchel dual ofbM- Equivalently, 
fu is given by: 

fM{x,x*) = sup{(a,x*) + (x,a*) - {a, a*) : {a, a*) G M} 

Proposition 5.2 (Properties of the Fitzpatrick function) Let M C X x X* be a 
graph. Then the associated Fitzpatrick function fM has the properties: 

(a) fM is convex and lower semicontinuous, 

(b) the graph M is maximal monotone if and only if: 

(bl) for any (x,x*) £ X x X* we have fM{x,x*) > {x,x*) and 
(b2) we have equality fM{x,x*) = {x,x*) if and only if {x,x*) G M. 



(c) if M is a BB-graph then the function 

fM{x,x*) , {x,x*) G doniM x imM 



1^ +00 otherwise 

is a strong bipotential. 

Proof. By construction any Fitzpatrick function is convex and lower semicontin- 
uous. For proving (b) it is enough to use the following characterization of maximal 
monotone graphs: M is maximal monotone if and only if 

(il) for any {x,x*) e X x X* we have 

mi{{x-a,x* -a*) : {a,a*) G M} < 

(i2) (x, X*) G M if and only if 

inf{(a;-a,a;*-a*) : {a,a*) G M} = 

The Fitzpatrick function associated to the graph M can be written as: 

f]\f{x,x*) = {x,x*) — ini{{x — a,x* — a*) : {a,a*) £ M} 

Therefore (bl), (b2) follow from (il), (12) respectively. 

In order to prove (c) we have to check (BIS), (B2S) definition 13.31 Let x G 
domM. Then 

m{{gM{x,x*) — {x,x*) : x* G imM} = 

= -inf{(x-a,a;* -a*) : {a, a*) e M , x* e imM} = 

The proof of (B2S) is similar. D 

Further we describe selfdual lagrangians, with the notations from Ghoussoub 
|18] . See the mentioned paper and the references therein, especially [T9] [20] for 
more informations on the variational theory associated to selfdual lagrangians. Here 
we point out that selfdual lagrangians are particular cases of bipotentials. 

Definition 5.3 Let X be a reflexive Banach space. To any function L G T(){XxX*) 
we associate the following operators: 

(i) 5L:X ^2^' defined by: 

5L{x) = {x* G X* : {x,x*) G dL{x,x*)} 

(here dL{x,x*) is the subdifferential of L), 

(ii) dL:X^ 2^* defined by: 

dL{x) = {x* eX* : L{x,x* = {x,x*)} 
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Definition 5.4 A function L e Tq{X x X*) is a selfdual lagrangian if L = L* . 

The following is a slight reformulation of lemma 2.1 and proposition 2.1 [18] 

Proposition 5.5 If L is a selfdual lagrangian such that for some xq £ X the func- 
tion L(xo, •) 'i'S bounded on the balls of X* , then L is a strong bipotential and we 
have M{L) = G{dL) = G{5L). Moreover, in this case for any x* G X* there exists 
X & X such that x* £ 6L{x) and 

L{x,x*) - {x,x*) = mi{L{x,x*) - {x,x*) : X e X} = 

Proof. As mentioned before, from lemma 2.1 [18j we get that for any selfdual 
lagrangian L we have dL = SL. By definition any selfdual lagrangian is convex, 
lower semicontinuous and for any (x,x*) G X x X* we have L(x,x*) > {x,x*), as 
a consequence of the Fenchel inequality in X x X*. The fact that L is a strong 
bipotential (conditions (BSl), (BS2) definition 13. 3p . as well as the final part of the 
conclusion are straightforward reformulations of the conclusion of Proposition 2.1 

The following proposition is an application of lemma 3.1 and proposition 3.1 [18]. 

Proposition 5.6 Let M C X x X* be a maximal monotone graph. Then there 
exists a selfdual lagrangian L^j such that G{5Lm) = M . 

Proof. Let /a/ be the Fitzpatrick function of M. Then, according to lemma 
3.1 [E] the Fitzpatrick function /m satisfies the hypothesis of proposition 3.1 [18j- 
Therefore the selfdual lagrangian defined by 

f 1 1 1 

Lm{x,x*) = inf <^ -/Af(xi,xi) + -/|j(x2,X2) + -||xi -X2IP+ 

+ glkl -Xlf ■■ {X,X*) = -{Xl,xl) + -{X2,XI) 

mentioned in the proof of the proposition 3.1 [T8] as "the proximal average" between 
fM and fl,^ is the one needed. D 

6 Separable bipotentials 

If (/> : X — > M is a convex, lower semicontinuous potential, consider the multivalued 
operator d(f) (the subdifferential of </>). The graph of this operator is the set 

M{ct)) = {{x,y)eXxY I 0(x) + r(y) = (x,y)} . (6.0.1) 

M{(j)) is maximal cyclically monotone [32] Theorem 24.8. Conversely, if M is closed 
and maximally cyclically monotone then there is a convex, lower semicontinuous (j) 
such that M = M{4>). 

9 



To the function (p we associate the separable bipotential 

b{x,y) =(^(x) + (l)*{y). 

Indeed, the Fenchel inequahty can be reformulated by saying that the function b, 
previously defined, is a bipotential. More precisely, the point (b) (resp. (c)) in the 
definition of a bipotential corresponds to (i) (resp. (ii)) from Fenchel inequality. 
The bipotential b and the function cj) have the same graph, that is M(b) = M{(p). 

7 Bipotentials for monotone, non maximal graphs 

The following two results are from the paper pT] (theorem 3.1 and corollary 3.3). 

In the theorem below it is shown that intersections of two maximal monotone 
graphs are sometimes representable by a bipotential. Therefore there exist bipoten- 
tials b with M{b) monotone, but not maximal. 

Theorem 7.1 Let bi and 62 be separable bipotentials associated respectively to the 
convex and lower semicontinuous functions (pi, (p2 : X ^>-MU {+00}, that is 

bi{x,y) = (piix) + (t)*{y) 

for any i = 1,2 and {x,y) G X x Y . Consider the following assertions: 

(i) b = max{bi, 62) is a strong bipotential and M{b) = M{bi) n M{b2)- 

(ii') For any y £ douKp* D dom(p2 and for any A G [0, 1] we have 

(XcPi + (l-A)<A2r(y) = XcPliy) + (1-A)0^(y) (7.0.1) 

(ii") For any x G dam (pi f] dom(p2 and for any A G [0, 1] we have 

(A0* + {l-\)(Piy{x) = A0i(x) + (l-A)</>2(x) (7.0.2) 

Then the point (i) is equivalent with the conjunction of (ii'), (H"), (for short: (i) 
■^^ ( (ii') AND (ii") ) ). 

In the proof of this theorem we make use of a minimax result by Sion |i43j . Notice 
that in section [T2t we use another minimax result of Fan [29] in the proof of theorem 
1221 

The conditions (ii'), (ii") from Theorem 17.11 imply relations which can be ex- 
pressed with the help of inf convolutions. Consider (pi,(p2 G ro(X). For any 
A G (0, 1) we introduce two functions defined on X by: 

fl,x{x) = Xm\^) , /2,a(x) = (1-A)02(: ^ ^ 



'A ' ' ■'^""^ ' ' '^^'l-X 
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Proposition 7.2 Let 01,02 S ^oi^) such that 

b{x,y) = max(0i(x) + 0i(y),02(x) + 02(y)} 

is a strong bipotential. Then, with the previous notations, for any x £ dam (pi f] 
dom(f)2 and for any A G (0, 1) the subdifferential of the inf- convolution /i, An/2, A has 
the expression: 

5(/i,An/2,A)(x) = dMx) n dM^) 



8 Bipotentials and inequalities 

A source of interesting bipotentials is provided by inequalities. Here we discuss 
about the Cauchy-Bunyakovsky-Schwarz inequality and about an inequality of Fan 
concerning eigenvalues of symmetric matrices. 

8.1 Cauchy bipotential 

Let X = y be a Hilbert space and let the duality product be equal to the scalar 
product. Then we define the Cauchy bipotential by the formula 

b{x,y) = \\x\\ \\y\\. 

Let us check the Definition (j3.ip The point (a) is obviously satisfied. The point 
(b) is true by the Cauchy-Bunyakovsky-Schwarz inequality. We have equality in the 
Cauchy-Bunyakovsky-Schwarz inequality b{x, y) = {x, y) if and only if there is A > 
such that y = Ax or one of x and y vanishes. This is exactly the statement from the 
point (c), for the function b under study. 

The graph of the Cauchy bipotential is not monotone. 

8.2 Hill bipotential 

Let S{n) be the space of n x n real symmetric matrices. There is a bipotential 
expressing that two matrices X and Y have a simultaneous ordered spectral decom- 
position, [44j . In Mechanics, a constitutive law between two tensors implying that 
they admit the same eigenvectors is said to be coaxial [39], j44j . 
The space S{n) is endowed with the scalar product 

{X,Y) = tr{XY) 

Consider the function b : S{n) x S{n) -^ M with the expression 

biX,Y) = Ai(X)Ai(y) + ... + An(X)A„(y) 
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where for any X G S{n) Ai(X) > X2{X) > ... > Xn{X) are the eigenvalues of X 
ordered from the largest to the smallest. With the help of this function we can write 
one of Fan's inequalities [28] as: for any X,Y £ S{n) we have 

b{X,Y)>{X,Y) 

with equality if and only if X and Y have the same eigenvectors with preservation 
of the order of the eigenvalues. 

In [44] is proved that for n = 3 the function 6 is a bipotential, called Hill 
bipotential due to applications dealt with by Hill in mechanics. A similar proof, 
involving majorization, can be done for the case of a general n. The graph of the 
Hill bipotential is not monotone. 

9 Bipotentials in non smooth mechanics 

A simple example of a monotone operator in non smooth mechanics is provided 
by the following model of plasticity of metals. X = Y is the space of n x n real 
symmetric traceless matrices with the pairing {x, y) = tr xy and the associated norm 
||x|| =1 {x, y) 1 2 . Let c be a non negative constant. The plasticity operator is defined 
by 

Tp{0)=K = {yeY\ \\y\\ <c} 

otherwise 

Tp (x) = f--7 
\\x\\ 

In plasticity, x is the plastic strain rate tensor, y is the deviatoric stress tensor 
and c is the yield stress. The closed convex set K is the plastic domain and the 
irreversible or plastic deformations varies when x ^ 0. The plastic model is not 
limited to the metals but can be used also for soil materials. 

A more involved example is the associated Driicker-Prager model where the 
variable x is the plastic strain rate tensor and y is the stress tensor (both seen as 
elements of S{n)). With usual notations, the tensors x and y are split into their 
deviatoric and spheric parts: 

x = x, + lx,I, y = y, + y,I 

where I is the identity operator, Xh = tr (x) and yh = ^tr (y). As the decomposition 
is unique, we can write in short x = (xd^Xh), y = {yd,yh) and the duality pairing is 
{x,y) =tr{xdyd) + Xhyh- 

The convex cone parameterized by the friction angle 99 G (Oif)! the cohesion 

stress c > 0, and r = 3\/2/V9 + 12tan^ if ([7j), with the vertex aX v = /, 

tan 93 

given by 

K = {y eY \ I \yd\ \<r{c- yh tan 99)} 
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is called the plastic domain. 

The multivalued operator corresponding to the associated Driicker-Prager model 
is defined by: Tjjp (0) = K, if Xh = rWxdW tamp then 

Top (x) = \y eY \3r]>0, y = v + ■>] i j — r^ v 

[ VlFdII re 

if Xh > '''lla^dll tamp then Tdp (x) = {v}, otherwise Tdp (x) = 0. 

9.1 Driicker-Prager non associated plasticity 

The non associated Driicker-Prager model is characterized, as previously, by 
the friction angle ip S (O, ^) , the cohesion stress c > 0, and r = 3-\/2/ y 9 + 12 tan^ ip 
but also by a new parameter, the dilatancy angle 6 G [0, (p). Once again, X = Y = 
S{n) and we use the splitting into deviatoric and spheric parts. The associated 
multivalued operator is Tna defined by: Tna (0) = K, if x^ = r\\xii\\ tan0 then 

Tna{x) = <y eY \3r]>0, y = v + r]i - — v 

I VlFrfll re 

if Xh > r\\xd\\ tanO then Tna {x) = {v} otherwise Tna (x) = 0. 

If we put 9 = ip then we recover the operator T^p of the associated Driicker- 
Prager model defined above. 

The non-associated Driicker-Prager law y £ Tna{x) is equivalent with the follow- 
ing differential inclusion: 

X + -{xh + r ||x|| (tan(/> - tan 6*)) G dxKiv) 

According to [9] [39] [25], this inclusion can be written as bp{x,y) = {x,y), for the 
bipotential 

bpi^^y) = XK{y) + XKp{x) + -^^^^^h + 

+ r \\x\\ (tan (f) — tan9)[yh — 

\ tan( 

The last term in this expression is a coupling term which gives the implicit character 
to the constitutive law. 

9.2 Coulomb's dry friction law 

Another interesting operator comes in relation with the unilateral contact with dry 
friction or Coulomb's friction model. T Consider two bodies f^i and ^2 which 
are in contact at a point M, with n the unit vector normal to the common tangent 
plane and directed towards Qi. The space X = M^ is the one of relative velocities 
between points of contact of two bodies, and the space Y, identified also to M'^, is 
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the one of the contact reaction stresses. The duahty product is the usual scalar 
product. We put 

{xn,xt)eX = RxR\ {yr„yt)eY = RxR\ 

where Xn is the gap velocity, xt is the sliding velocity, y„ is the contact pressure and 
yt is the friction stress. The friction coefficient is // > 0. 

The graph of the law of unilateral contact with Coulomb's dry friction is the 
union of three sets, respectively corresponding to the 'body separation', the 'sticking' 
and the 'sliding'. 

M = {{x,0) e X xY I x„<0}U (9.2.1) 

U {{0,y)£XxY I ||yt ||<//y„}U 
U <{x,y) G X xY I Xn = 0, xt^O, yt = fJ-yni 



It is well known that this graph is not monotone, then not cyclically monotone. 
This law can be written as the following differential inclusion ([33] [3l] [IQ] |39j): 

{xn - n\\xt\\)n + xt e dxK^iy) 

Let us consider the conjugate cone of the Coulomb's cone: 

K*^ = {{Xn,Xt) £ X I /U II xj II +x„ < 0} . 

We shall use also a second pair of conjugate cones: 

Ko = {(yn,0) G y I yn > 0} , K^ = {(x„,Xi) eX \ Xn < 0} . 

The graph of the law of unilateral contact with Coulomb's dry friction is the 
graph of the following bipotential, previously given in [33j : 

b{x,y) = Wn II Xt II +XKf,{y) + xk^,{x) ■ 

9.3 Coaxial laws 

Consider the non monotone operator defined by Ti^o (0) = M", otherwise 

T^so (x) = {y G y I 3A > 0, y = Ax} 

An operator S : M" -^2 is strongly isotropic if its graph is contained in the graph 
of TisQ. The graph of the non monotone operator Tiso is the graph of the Cauchy 
bipotential. 

The eigenvalues of any matrix x E S{n) can be conventionally ordered from the 
largest to the smallest: Ai(x) > A2(x) > ... > A„(x). Two elements x,y £ S{n) 
are said coaxial if they have the same eigenvectors with preservation of the order 
of the eigenvalues: the largest eigenvalue of x is associated to the eigenvector of y 
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corresponding to the largest eigenvalue of y, and so on. An operator T : S{n) -^ 2^^"'' 
is coaxial ([I2],[ll]) if for any y £ T{x), y and x are coaxial. The graph of any coaxial 
operator is contained in the non monotone operator defined by Th (0) = S{n), 
otherwise 

Th (x) = {y GY \ X and y are coaxial} 

We have seen in section 18.21 that the graph of a coaxial operator is included in the 
graph of the Hill bipotential. 

10 Numerical methods based on the bipotential frame- 
work 

In applications bipotentials are interesting because of the associated implicit nor- 
mality rules. The properties of bipotentials allow to discretize an evolution problem 
into a series of minimization problems concerning the associated bifunctional. 

For example let us consider a body with reference configuration 0, with the 
boundary decomposition 9$7 = Fi U r2 U Fs. At a fixed moment t we have imposed 
velocities iit on the part Fi of its boundary, imposed forces ft on F2 and on F3 the 
body is in unilateral contact with friction with a rigid foundation. The unit normal, 
pointing outwards, of the boundary dfl is denoted by n and for any stress field a in 
ri we denote an = crn. 

Suppose that the body is made by a plastic material described by a bipotential 
bp dependent on the strain rate e(n) and the stress a. The contact with friction is 
described by the bipotential be- associated with the Coulomb law. 

For any pair (v, r) of kinematically admissible velocity field v and statically 
admissible stress field r we introduce the bifunctional 

B{v,t) = / bp{e{v),T) dx + / bc{-v,Tn) ds - 
JQ Jr-j, 

- / Tn-Ut ds - / ffV ds 

Then, by using integration by parts and the properties of the bipotentials which are 
involved, one can show that 

B{v,t) > 

and that B{u, o") = if (n, a) is the pair formed by the velocity field and associated 
stress field of the body at moment t. 

Therefore we may try to numerically minimize the bifunctional in order to find 
the solution of the (quasistatic) evolution problem. Alternatively, the bifunctional 
can be adapted in order to use a predictor/corrector scheme. 

As applications we may cite the bound theorems of the limit analysis ( |38j , [8] ) 
and the plastic shakedown theory ([^, [12], [9j, [6j) which can be reformulated by 
means of weak normality rules. The bipotential method suggests new algorithms, 



15 



fast and robust, as well as variational error estimators assessing the accurateness of 
the finite element mesh ([22], [23], [37], [ID], [7], [2S], [20]). Such algorithms have 
been proposed and used in applications to the contact mechanics |13j . the dynamics 
of granular materials ([15]. [16j . [IZ][12]), the cyclic plasticity of metals [37] and the 
plasticity of soils ([^, [25]). 

A very challenging subject seems to be the extension of the mathematical results 
of Ghoussoub [19] [20] for the particular case of selfdual lagrangians, to general 
bipotentials, or at least to bipotentials constructed from bipotential convex covers. 

11 Construction of bipotentials 

Let Bp{X, Y) be the set of all bipotentials b : X xY ^ M.U {+00} . We shall need 
the following notion of implicitly convex functions. 

Definition 11.1 Let A be an arbitrary non empty set and V a real vector space. 
The function f : A x V ^ M is implicitly convex if for any two elements 
(Ai, zi), (A2, 2:2) G A X y and for any two numbers a, fi £ [0,1] with a + (3 = 1 
there exists A E A such that 

f{X,azi+ (3z2) < af{Xi,zi) + (3f{X2,Z2) ■ 



In the following bipotential convex covers are defined, as in definition 4.2 

Definition 11.2 A bipotential convex cover of the non empty set M is a func- 
tion A G A I— > 6>, from A with values in the set Bp{X, Y), with the properties: 

(a) The set A is a non empty compact topological space, 

(b) Let /:AxXxy^RU {+cxd} be the function defined by 

f{\x,y) = bx{x,y). 

Then for any x £ X and for any y £ Y the functions f{-,x, •) : A x y ^ M 
and f{-,-,y) : A X A ^ M are lower semi continuous on the product spaces 
A X y and respectively A x A endowed with the standard topology, 

(c) We have M = \J M{bx). 

(d) the functions f{-,x, ■) and /(•, -jy) are implicitly convex in the sense of Defi- 
nition \lin\ 
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This notion generalizes the one of a bi-implicitly convex lagrangian cover. 

see Definitions 4.1 and 6.6 [TU]. Here we shall give only the definition of a convex 
lagrangian cover, without the bi-implicit convexity hypothesis. 

Definition 11.3 Let M C X x Y be a non empty set. A convex lagrangian 

cover of M is a function X £ A ^^ (px from, A with values in the set of lower 
semicontinuous and convex functions on X, with the properties: 

(a) The set A is a non empty compact topological space, 

(b) Let f: AxXxY^M be the function defined by 

f{\x,y) = (t>x{x) + (t>l{y). 

Then for any x £ X and for any y GY the functions f{-,x, •) and /(•, •, y) are 
lower semicontinuous from A with values in the set of lower semicontinuous 
and convex functions on X, endowed with pointwise convergence topology, 

(c) we have 

M = [j M{<Px) . 
AeA 

A bipotential convex cover A G A i-^ 6a such that for any A G A the bipotential 
b\ is separable is a bi-implicitly convex lagrangian cover, see Definitions 4.1 and 6.6 
|10j . . For such covers the sets M[b\)&Te maximal cyclically monotone for any 
AeA. 

General bipotential convex covers are not lagrangian. (see remark 6.1 [10] for 
a justification of the "lagrangian" term). In the language of convex analysis this 
means that in general the sets M{b\) are not cyclically monotone. An example is 
given in section 5 [11], of a bipotential convex cover of the graph of the Coulomb's 
dry friction law, which is made by monotone, but not maximally monotone graphs. 

In sections 5 and 8 [10] it is explained that not any BB-graph admits a convex 
lagrangian cover. Moreover, there are BB-graphs admitting (up to reparametriza- 
tion) only one convex lagrangian cover, as well as BB-graphs which have infinitely 
many lagrangian covers. The problem of describing the set of all convex lagrangian 
covers of a BB-graph seems to be difhcult. 

Definition 11.4 Let X ^^ bx be a bipotential convex cover of the BB-graph M . To 
the cover we associate the function b : X x Y ^ M.U {-|-cxd} by the formula 

b{x,y) = ini {bx{x,y) .■ A G A} 



We obtained in theorem 4.6 [TT] the following result. 
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Theorem 11.5 Let X >-^ (px be a bipotential convex cover of the BB- graph M and 
b : X X Y ^ R the associated function introduced in Definition \ll-4\ Then b is a 
bipotential and M = M(b). 

In the case of M = M{(j)), with (j) convex and lower semi continuous (this cor- 
responds to separable bipotentials), the set A has only one element A = {A} and 
we have only one potential (p. The associated bipotential from Definition 111.41 is 
obviously 

b{x,y) = cpix) + (t)*{y) . 

This is a bipotential convex cover in a trivial way; the implicit convexity conditions 
are equivalent with the convexity of 0, (p* respectively. 

12 One more construction result 

For simplicity, in this section we shall work only with lower semicontinuous convex 
functions cf) with the property that (p G r(X, M) and its Fenchel dual (j)* e r(y, M). 
We reproduce here the following definition of convexity (in a generalized sense) , 
given by K. Fan [29] p. 42. 

Definition 12.1 Let X, Y be two arbitrary non empty sets. The function f : 
XxY ^fM. is convex on X in the sense of Fan if for any two elements xi, X2 ^ X 
and for any two numbers a, (3 £ [0, 1] with a + /3 = 1 there exists a x £ X such that 
for all y £Y : 

f{x,y) < af{xi,y)+Pf{x2,y). 



With the help of the previous definition we introduce a new convexity condition 
for a convex lagrangian cover. 

Definition 12.2 Let X >-^ (px be a convex lagrangian cover of the BB-graph M , in 
the sense of definition \11.3[ Consider the functions: 

g-.XxAxX^m , h:YxAxY^R , 

given by g{x,X,z) = cpxix) -(j)x{z), respectively h{y,X,u) = (pKy) - (pKu). 

The cover is Fan bi-implicitly convex if for any x £ X , y gY , the functions 
g{x, ■, •), /i(y, •, •) are convex in the sense of Fan on A x X , A x Y respectively. 

Recall the following minimax theorem of Fan [29], Theorem 2. In the formulation 
of the theorem words " convex" and " concave" have the meaning given in definition 
112.11 (more precisely / is concave if — / is convex in the sense of the before mentioned 
definition) . 
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Theorem 12.3 (Fan) Let X he a compact Hausdorff space and Y an arbitrary set. 
Let f be a real valued function on X xY such that, for every y £ Y, /(•, y) is lower 
semicontinuous on X. If f is convex on X and concave on Y , then we have 

min sup /(rr, y) = sup min /(x, y) . 
xGX y£Y y&y ^'S-^ 



The difficulty of theorem ll 1.51 boils down to the fact the class of convex functions 
is not closed with respect to the inf operator. Nevertheless, by using Fan theorem 
I12.3l we get the following general result. 

Theorem 12.4 Let K he a compact Hausdorff space and X >-^ (j)x £ T{X, M) be a 
convex lagrangian cover of the BB-graph M, with (j)\ € r(y, M) for any A G A, such 
that: 

(a) for any x £ X and for any y £ Y the functions A B X ^^ 4>x{x) G I^ o-i^-d 
A 3 A 1-^ (/>^(y) G M are continuous, 

(h) the cover is Fan hi-implicitly convex in the sense of definition \12.2l 

Then the function b : X x Y ^M. defined by 

b{x,y) = inf {(t>x{x) + ^liy) | A G A} 

is a bipotential and M = M{h). 

Proof. For some of the details of the proof we refer to the proof of theorem 4.12 
lll.5l in [lOj . There are five steps in that proof. In order to prove our theorem we have 
only to modify the first two steps: we want to show that for any x £ dom{M) and 
any y £ im(M) the functions b{-,y) and b(x, •) are convex and lower semi continuous. 
For {x,y) £ X xY let us define the function xy : A x X ^ M by 

xy{X,z) = {z,y) + (l)x{x) -(px{z) . 

We check now that xy verifies the hypothesis of theorem ll2.31 Indeed, the hypothesis 
(a) implies that for any z £ X the function xy{-, z) is continuous. Notice that 

xy{X, z) = (z, y) + g{x, A, z) . 

It follows from hypothesis (b) that the function xy is convex on A in the sense of 
Fan. 

In order to prove the concavity of 'xy on X, it suffices to show that for any 
zi,Z2 £ X, for any a, /? G [0, 1] such that a + /3 = 1, we have the inequality 

xy(A, azi + I3z2) < axy{X, zi) + pxy{X, Z2) 
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for any A G A. This inequality is equivalent with 

{azi + f3z2,y) - 4)x{azi + (3z2) < a{{zi,y) - (pxizi) + f3{{z2,y) - (t)\{z2) 

for any A G A. But this is implied by the convexity of (px for any A G A. 

In conclusion the function xy satisfies the hypothesis of theorem ll2.31 We deduce 
that 

minsupxy(A, z) = supminxy(A, z) . (12.0.1) 

Let us compute the two sides of this equality. 

For the left hand side term of (jl2.0.ip we have: 

min sup xy(A, z) = min sup {{z, y) + (px{x) — (pxi^)} = 
AgA ^gx -^gA 2gx 

= min<^ (pxix) + sup{(2;,y) - (pxiz)} 
agA [ 2gx 

= min {(f) x{x) + (j)*x{y)} = b{x,y) . 
AeA 

For the right hand side term of (jl2.0.ip we have: 

sup minxy(A, z) = sup min {(z, y) + (pxi^) — 4'x{z)} 
z€X AeA zex aga 



= sup <^ {z, y) - max {(pxiz) - (px{x)} 

zeX [ A6A 

Let rr : X ^ M be the function 

x{z) = max{cl)x{z) - (pxix)} . 
AeA 

Then the right hand side term of (|12.0.ip is in fact: 

supmin^(A,z) = x*(y) . 

zex AeA 

Therefore we proved the equality: 

b{x,y)=x*{y) . 

This shows that the function b is convex and lower semicontinuous in the second 
argument. 

In order to prove that b is convex and lower semicontinuous in the first argument, 
replace (px by (f>*x in the previous reasoning. D 
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